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Abstract. Lorentz's reciprocity lemma and Feld-Tai reciprocity theorem show 
the effect of interchanging the action and reaction in Maxwell's equations. We 
derive a free-space version of these reciprocity relations which generalizes the 
conservation of the momentum-energy tensor. This relation corresponds to 
the interference conservation of electromagnetic waves. We show that for any 
transformation or symmetry that leaves Maxwell's equations invariant, we can 
modify the reciprocity relation to introduce a conserving density, optical flux and 
stress tensor extending Noether's theorem to a different context. We apply this 
method to transformations that can be expressed as Hermitian operators and 
more specifically, we define the operators associated with the optical energy, spin, 
linear and angular momentum. 
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1. Introduction 

Defining the optical spin and angular momentum density of the electromagnetic 
field is an important aspect in the study of optical vortices and their evolution and 
dynamics. Indeed, optical vortices are associated with a circular flow of energy, but 
how can we deflne a physically meaningful quantity that helps us understand this 
flow? One possible approach is to use operators taken from Quantum Mechanics on 
the electromagnetic fields. This can be done within the paraxial approximation pn,'2] or 
based on the second quantization [Sj. Within this paraxial approximation, Maxwell's 
equations become scalar and are formally identical to Schrodinger's equations making 
it possible to use quantum mechanical operators. Unfortunately, defining such 
quantities for vectorial fields is not as straight forward. Various similar methods 
have been proposed in the literature based on the angular momentum of the Poynting 
vector [H H]. More recently, the energy fiow around vortices [6j has been equated 
to optical currents by invoking the "democratic" principle between the electric and 
magnetic fields [7]. It is this last approach that we build upon and generalise in 
this paper by showing crucially that this optical current is only one part of whole 
conserving energy-momentum tensor. The method presented here is based on the 
electromagnetic reciprocity theorem and the superposition principle. The generality 
of this approach makes it applicable to other conserving currents, each corresponding 
to a symmetry, and even to other field evolution equations. Finally, we conclude by 
discussing the link between these conserving currents, the symmetry operators and 
the quantum mechanical operators. 

2. Free space reciprocity relations 

The starting point of our paper is the optical reciprocity theorem, which defines 
Newton's third law for optical fields [3[9]. More generally, the reciprocity relations 
describe what happens when cause and effect are exchanged. For Maxwell's equations, 
the reciprocity theorem takes two forms. The Lorentz reciprocity lemma which links 
the electric field created by an emitting current and the electric field from the receiving 
current and the Feld-Tai reciprocity theorem which links the same currents but via 
the generated magnetic fields |10[ lllj . Here, we are not interested in the current- 
field interaction part of the reciprocity relation but only in the electromagnetic free- 
field behaviour within the reciprocity theorem. In this case, the reciprocity theorems 
simplify to an interference conserving relation linking two independent solution of 
Maxwell's equations. In the appendix, we derive this reciprocity relation in the case 
of the Riemann-Silberstein form of Maxwell's equations [T2t[T3] that formally is similar 
to Schrodinger's equations. 

Let us consider Maxwell's equations in free space, 

V • eoE = 0, 

V • MoH = 0, 

V X E = -fiodtil, ^ ' 

V X H = eo^tE, 

where E and H are the complex electric and magnetic vector fields and where eo and no 
the vacuum permittivity and permeability. The energy and momentum associated with 
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these fields is given by E the energy density, S the Poynting vector and a Maxwell's 
stress tensor [14] 

E{J^) ^i(eoE*.E + ^oH*-H), 

Sm ==i(E* X H + E X H*), 

2^ (2) 



= — (^(eoE* • E + ^oH* • H) / - eoE* ® E 
-eoE ® E* - ^oH* ® H - ^qH (g) H*) , 

where (g) stands for the tensor product, and * stands for the complex conjugate, 
c = l/y/iojM) for the speed of light, / for the identity 3x3 matrix and T for the pair of 
vectors ~ (E, H). The conservation of energy and momentum of the electromagnetic 
field is given by: 

V.S(^) + 5*£;(^) = 0, 

V-CT(^) + 5tS(^) = 0. ^ ' 

An intuitive way to deduce the free space reciprocity relations is through the 
use of Hnear superposition [9]. Here, we have simplified this approach by considering 
only the linear superposition T = T\ + T2 of two fields T\ and T2, both solutions of 
Maxwell's equations H]). As the energy, momentum and stress tensor are all quadratic 
in the field, each of them can be split these into three parts. Two parts involve only 
terms in T\ and Ti respectively while the remaining terms involve cross products 
between these fields. This decomposition can the represented as: 

S(.F) =S(.Fl) + Si2(.Fi,.F2) + S(.F2), (4) 

where the expression with the subscript 12 groups all the terms involving cross products 
between the two fields. These expression are given by: 

^12(^1,^2) (eoEt • E2 + /^oHt • H2) , 
Si2(.Fi,.F2)=i(EI xH2 + E2XHt), 

(5) 



^12(^1,^2) - y ((eoE^ • E2 + MoH^ • H2) / - eoEJ ® E2 
-eoE2 «) E^ - [IqI^I (g) H2 - /ioH2 ® H^) . 

These expressions describe the energy and momentum associated with the interference 
of the two fields and correspond to its conservation. Indeed, each of the two fields T\ 
and Ti fulfill the conservation relations ([3]). Therefore, when substituting Q into ^ 
we find: 

V-Si2(.Fl,.F2)+9tSl2(.Fi,.F2) =0, 

V-S?i2(.Fl,.F2)+9tSi2(.Fl,.F2) =0. ^ ' 

These relations describe the conservation of the interference of two solution of 
Maxwell's equations and represent the free-space equivalent to Lorentz lemma [8] 
and the Feld-Tai reciprocity relations |10[ lllj . Physically, relation ([6]) means that 
if the superposition of two fields is associated with a specific interference energy 
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density, then energy density this will be conserved regardless of the linear propagation 
transformation acting on each of the fields. Otherwise stated, it is not possible to 
treat the two field distributions independently of each other which is equivalent to 
a classical energy entanglement (property that has also been observed in |T5|). In 
the appendix, we present the same procedure for the Riemann-Silberstein form of 
Maxwell's equations. 

3. Invariant transformations and symmetries 

In this section, we consider the effect of a field or space transformation, which leaves 
Maxwell's equations invariant, on the free space reciprocity relation defined in the 
second section. These transformations are each Hnked to a symmetry of Maxwell's 
equations and we show that for each such symmetry there exists a conserving optical 
current or more generally a canonical conserving energy-momentum tensor. Our result 
extends Noether's theorem [16] to the reciprocity theorem offering a direct way to 
deduce symmetric canonical energy-mometum tensors from each other (see [17] for a 
discussion of the symmetry problem). This first result of the paper is stated in the 
form of a theorem: 

If J- is a solution of Maxwell's equations, T a transformation that leaves Maxwell's 
equations invariant andTT is the transformed field still solution of Maxwell's equation, 
then Ei2 {T,TT) , S 12 (.?-", and ai2 {T,TT) define the canonical energy density, 
current and stress tensor associated with this transformation or symmetry. These 
quantities form the canonical energy-momentum tensor and verify the conservation 
relation 

This theorem can be proved directly by substituting T\= T and T2 — TJF in the 
reciprocity relations |[5|) and This means that for every invariant transformation 
of Maxwell's equations, we can construct a conserving quantity together with its 
conserving current. This result is related to Noether's theorem that uses invariant 
symmetries of the Lagrangian to deduce, for each symmetry, a canonically conserving 
current [16j . 

Additionally for each symmetry or invariant transformation, we can define the 
following primary integral of the system: 



These conserving canonical quantities can be associated with each solution and 
symmetry of Maxwell's equation and are constant in time. Here, we note that these 
quantities are not linear with the respect to the field but quadratic. A prime example 
of such a conserving quantity is the energy of the electromagnetic wave which is 
< Ei2{J^) >i where the transformation corresponds to the identity transformation I. 

4. Canonical superposition principle 

In the fourth section, we deal with the interplay between conservation, superposition 
and interference of independent solutions of Maxwell's equations. Indeed, a 
superposition of multiple such solutions is still a solution. However, because of 
the interference between these different solutions, the conserving quantities are not 




(7) 
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necessarily additive. Starting from this observation, we formulate the principle 
of canonical superposition implying a simultaneous linearity in the fields and in 
the conserving quantities in a Hnear superposition. To satisfy this principle, the 
electromagnetic field needs to be decomposed onto the eigenfunctions of the operators 
associated with each symmetry. Additionally, for Hermitian operators we show that a 
complete sets of pair-wise commuting operators correspond to compatible conserving 
quantities. However, there are no electromagnetic fields that simultaneously satisfy 
the canonical superposition principle for non-commuting operators. This approach 
enables the introduction of the photon wave function fL8\ and a first quantisation of 
Maxwell's equation |19j . 

In this section, we restrict ourselves to invariant transformations that are linear 
and can be expressed as a linear operator A acting on the complex vector field 
defined earlier. The second result of this paper can be written in the form of the 
principle: 

For a vector field to be observed as an independent entity when interfering with 
another external field, its conserving canonical quantities must be additive to those of 
the external field. 

Mathematically this can be expressed in the following way. Let the total field 
be decomposed into a superposition of fields = ^ . where each !Fi is a solution of 
Maxwell's equations. For each part of the field J^i to be observed as an independent 
entity, the total conserving canonical quantity must be a sum of all the individual 
conserving canonical quantities 



Considering linear Hermitian operators A that leave Maxwell's equations 
invariant, we can verify that their eigenfunctions are fulfilling the canonical 
superposition principle. Indeed, we have: 



where Xi is the eigenvalue of XiTi — ATi . This is due to the fact that the eigenfunctions 
of a Hermitian operator are orthogonal and as such their interference integrated over 
the whole space cancels out. 

To appreciate that the canonical superposition principle is not automatically 
satisfied, we can look at the example of adding two identical electromagnetic fields 
similar to the process of stimulated emission. This superposition increases the total 
energy by a factor of 4. For the canonical superposition principle to be fulfilled with 
two identical fields, the two fields need to be in phase quadrature with respect to each 
other. In this case, adding the two fields together also adds their total energy. This 
means that the field associated with the stimulated emission needs to be out of phase 
with respect to the incident field such that its energy and field distribution are adding 
correctly. 

Finally, considering the simultaneous canonical superposition principle of multiple 
conserving quantities, we remark that there are quantities that are mutually exclusive. 
Indeed, for the superposition to work simultaneously for multiple Hnear operators, the 
eigenfunctions of these operators need to be shared. This is only possible for an 
ensemble of Hermitian operators that commute pair-wise. The maximum number of 




(8) 



< Ei2{:Fi) >A . 



< Ei2{T) >A 



^A, < Ei2{:Fi) >A 



(9) 
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quantities that fulfil the canonical superposition principle is given by a complete set 
of commuting operators. This also implies that there are conserving quantities that 
are mutually exclusive i.e. it is not possible to define a solution of Maxwell's equation 
that, at the same time, can be treated as an entity, conserving in a superposition, 
these canonical quantities. 

5. Applications 

In the final section, we apply the canonical superposition principle to define the spin, 
dynamic energy, linear and angular momentum currents. Four operators associated 
with these quantities form a complete set of pair-wise commuting operators with the 
non-diffracting Bessel beam as a common eigenfunction. We also verify the ratio 
between the dynamic energy and spin to be equal to the optical frequency [2Qj. The 
form of the spin current and density derived using this method is proportional to that 
found in the literature [H [7] . 

5.1. Classical energy-momentum 

Here, we remark that deriving the classical energy-momentum tensor from the free 
space reciprocity relations is straightforward. The transformation needed in this case 
is the identity transformation and implies the following relations: 

A = I, 

E,2 {T, AT) =i (eoE* • E + a^qH* • H) , 

Si2(.F,A.F)=i(E* xH + ExH*), (10) 

5i2 {T, AT) = y (^(eoE* • E + a^oH* • H) / - eqE* ® E 
-eoE «) E* - ^oH* (g) H - ^qH (g) H*) . 

5.2. Dynamic energy-momentum 

The first non-trivial Hermitian operator that leaves Maxwell's equations invariant is 
the derivative with respect to the time co-ordinate. In Quantum Field theory the 
time derivative operator corresponds to the invariance to the time translation and is 
associated with the energy of the system [16]. Using this derivative, all the static parts 
of the electromagnetic fields are eliminated and only the dynamic parts remains. The 
imaginary coefficient is necessary to maintain hermiticity of the dynamic operator and 
the energy-momentum tensor is given by: 

At = idt, 

Ei2 {T, AtT) =- (eoE* • 9tE -f y.^YL* ■ dtU) , 

Si2 (T, AtT) - ^ (E* X dtU + dt-E X H*) , (11) 

- 2 

5i2 (T, AtT) = ^ ((eoE* • S^E + ^ioH* • dtU) I- eoE* (g ^tE 
-eodtE ® E* - /ioH* ^ dtU - fiodtU «) H*) . 
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We remark that the eigenfunctions of the dynamic operator idt are 
monochromatic waves with an optical frequency defined as uj. For these waves, the 
dynamic energy-momentum tensor is proportional to the classical energy ifTOj) and to 
UJ. Unit-wise, the energy units can be recovered by making the electric and magnetic 
fields unit-less and using h as the conversion factor. Finally, using this conversion 
factor we can also introduce the momentum operator Ap = ihV which corresponds in 
Quantum Field theory to the translation independence. 



5.3. Optical spin flux 

The optical spin is associated with the circular polarisation of the light field. Here, we 
look for an Hermitian operator whose eigenfunctions are circular polarised light fields, 
which can be written as 

iZo 

-i/Zo 

where Zo = ^qc is the vacuum impedance. This operator corresponds to the duality 
transformation which leaves Maxwell's equations invariant. We also remark that the 
vacuum impedance does not break the Hermiticity of this operator as its purpose is to 
convert between the electric and magnetic field. The corresponding energy-momentum 
tensor is given by: 

^12 {T, ksT) =^ (E* • H - H* • E) , 

zc 

Si2 {T, k,T) = - (eoE* X E + AioH X H*) , ^^^^ 

CT12 {T, K,T) = - (^(E* • H - H* • E) / - E* (8) H 

-H «) E* + H* (8) E + E H*) . 
This describes the conservation of the circular polarisation as the light field 
propagates. Indeed, all the terms in the optical spin energy-momentum tensor are 
zero for linearly polarised light. Additionally, a circularly polarised plane wave 
has a constant optical spin density and the sign of this density depends on the 
polarisation handedness. The conservation of this definition of the optical spin has 
been theoretically asserted [H [U but not its associated stress tensor meaning it 
is not Lorentz invariant. Finally, the ratio between the optical spin density and the 
dynamic energy density as defined by (fTTj) is, as expected, u |20) . 



5.^. Optical angular momentum 

The angular momentum operator takes advantage of the cylindrical symmetry. Here, 
we have expressed it in cylindrical co-ordinates and we consider only a momentum in 
the z direction where (p is the azimuthal angle: 

= idcj,, 

{T, k^T) (eoE* • a^E + a^oH* • S^H) , 
Si2 {T, K^T) = ^ (E* X a^H + a^E x H*) , (13) 

ai2 = ^ ((eoE* • a^E + ^ioH* • ^^H) 7- cqE* ® ^^E 

-eoS^E (g) E* - /L(oH* ® d^H - fiod^^U «) H*) . 
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Higher order Laguerre-Gaussian and Bessel modes having a phase-front proportional 
to exp{il(j)) are eigenfunctions of this operator and their conserving canonical quantity 
corresponds to their azimuthal number I. Here again, making the electromagnetic field 
dimensionless and using Planck's constant h gives us the correct units for the angular 
momentum density. Additionally, the full momentum operator can be obtained by 
using Al = ihr x V which corresponds to a condensation of three operators that are 
non commutating between them. 

An interesting decomposition of the electromagnetic field can be achieved when 
considering the last three canonical conserving quantities. Indeed, the z-directional 
angular momentum operator together with the dynamic energy Aj, the z- 
momentum A^ and the spin Ag operators commute pair-wise and define a complete 
base of eigenfields. These vectorial electromagnetic fields correspond to circular 
polarized Bessel modes indexed using their azimuthal order, circular polarisation, 
longitudinal wave vector and optical frequency. 

6. Conclusion 

By using the free-space reciprocity relation, we have introduced a method to generate 
conserving relations for every invariant transformation of Maxwell's equations. For 
Hermitian operators, we have shown that their eigenfunctions decompose any 
electromagnetic field into a superposition, where each part can be distinguished by 
its canonical conserving quantities. The choice of these quantities is restricted by 
the number of Hermitian operators that commute pair-wise. Finally, we appHed this 
approach to define the dynamic energy momentum, the spin and angular momentum 
of the electromagnetic field. This results makes a direct first quantisation of Maxwell's 
equations possible and brings the concept of the classical electromagnetic photon wave 
function closer. Finally, by analogy it gives an insight into the real and imaginary part 
of the complex wave function in Quantum Mechanics (see Appendix) . 
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Appendix: Riemann-Silberstein equations and their reciprocity relation 

In vacuum. Maxwell's equations ^ are equivalent to: 



where the complex vector field is defined as F = (y^E + Zy^H)/2 and where electric 
and magnetic fields considered here are purely real. These are the Riemann-Silberstein 
form of Maxwell's equations[T2l [13] • We can determine the reciprocity theorem for 
these equations by considering the interference currents and density for a superposition 
of two fields F = Fi +F2, 



V • F = 







V X F = -dtF, 



(14) 



c 



Ers (Fi,F2) = F^ -Fs + Fi-F;, 

Srs (Fi, F2) = -ic{¥l X F2 -I- F* X Fi) 

GRs (Fi, F2) = -c^(F* (g) F2 -I- F2 «) F* H- Fi (g) F; -I- F; (g) Fi - EI), 



(15) 
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where Fi and F2 are each solutions of equation (fT4|) . Choosing Fi — F2 in these 
definitions we find the energy, the flow and stress tensor associated with equations lfT4|) . 

The two conservation relations of the interference currents become then the 
reciprocity relations for two solutions Fi and F2 of (fT4|) : 

V • Srs (Fi, F2) + dtEiis (Fi, F2) = 0, 

V-S?fls(Fi,F2) + atSfls(Fi,F2) = 0. ^ 
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